February 4, 2008 
BARI-TH/99-353 



Earth tides and Lense-Thirring effect 



Lorenzo lorio 



Dipartimento di Fisica dell' Universita di Bari, via Amendola 173, 70126, Bari, Italy 



Keywords: Tides, Node, Perigee, LAGEOS, LAGEOS II 



Abstract 

The general relativistic Lense-Thirring effect can be measured by inspecting a suitable combination of 
the orbital residuals of the nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS II. The solid 
and ocean Earth tides affect the recovery of the parameter by means of which the gravitomagnetic force 
is accounted for in the combined residuals. Thus an extensive analysis of the perturbations induced 
on these orbital elements by the solid and ocean Earth tides is carried out. It involves the / = 2 terms 
for the solid tides and the / = 2, 3, 4 terms for the ocean tides. The perigee of LAGEOS II turns out 
to be very sensitive to the I = 3 part of the ocean tidal spectrum, contrary to the nodes of LAGEOS 
and LAGEOS II. The uncertainty in the solid tidal perturbations, mainly due to the Love number k2, 
ranges from 0.4% to 1.5%, while the ocean tides are uncertain at 5% — 15% level. The obtained results 
are used in order to check in a preliminary way which tidal constituents the Lense-Thirring shift is 
sensitive to. In particular it is tested if the semisecular 18.6-year zonal tide really does not affect the 
combined residuals. It turns out that, if modeled at the level of accuracy worked out in the paper, 
the I = 2,4 m = and also, to a lesser extent, the I = 3 m = tidal perturbations cancel out . 
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1 Introduction 

If the Earth was perfectly spherical, according to Newton its gravitational potential would be 
GM®/r, where G is the Newtonian gravitational constant and is the mass of the Earth; 
the path of any artificial satellite in orbit around it would be a Keplerian ellipse with the Earth 
in one of its foci. This ellipse would remain fixed in inertial space and neither its shape nor 
its size would change; it could be parameterized by means of the so called Keplerian orbital 
elements [Sterne, 1960; Allison and Ashby, 1993]. In reality, the behavior of an object orbiting 
the Earth is much more complicated because our planet is not perfectly spherical due to several 
factors. First of all, the centrifugal force caused by its diurnal rotation makes the Earth oblate; 
its gravitational field, called geopotential, is no longer central but must be developed in a 
multipolar harmonic expansion \\Kaula, 1966|| : 



U = ^^{1 + E E (— )'P^n^(sin0)[a™cosmA + Si^smmX]} (1) 

' 1=2 m=Q ' 

where i?® is the Earth's mean equatorial radius, Pim is the associated Legendre function of 
degree / and order m, and and A are terrestrial latitude and east Greenwich longitude. Cim 
and Sim are the adimensional Stokes or geopotential coefficients [Lemoine et al, 1998] . Second, 
the shape of the Earth changes also in time due to the action of the solid, ocean and atmospheric 
tides which periodically lift the surface, tilt the vertical and redistribute the masses on our 
planet. This effect can be accounted for in a similar fashion to the static, centrifugal part of 
the geopotential given by eq.(|l]) introducing in it time varying-coefficients \\Eanes et al., 198^ 



The resulting effect of the static and dynamical non-sphericity of the Earth is that the path 
of any artificial satellite is more or less modified with respect to the unperturbed Keplerian 
ellipse. It is possible to account for these changes by assuming that the perturbed orbit can be 
still parameterized in terms of Keplerian elements which nevertheless slowly vary in time. This 
is so because the leading or point-mass term Uq = GM(^/r is predominant over the perturbing 
terms [Casotto, 1989]. Obviously, there is also a wide class of non gravitational effects [Milani 
et al., 1987] which contribute to modify the satellites' orbits; the concept of osculating ellipse 
fits also these cases and, in general, any situation in which a physical acceleration, small with 
respect to the Newtonian monopole, acts on the test body. 
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It should be pointed out that the experimental accuracy of techniques such as the Satellite 
Laser Ranging (SLR) has astonishingly grown up in the last decade: for example, today it is 
possible to measure the amplitude of a periodical perturbation on the node Q of LAGEOS with 
an accuracy of some milliarcseconds (mas in the following) . The accuracy for the same kind of 
measurement on the perigee cu is slightly worse (F. Vespe, private communication, 1999). This 
forces the researcher who wants to investigate some particular physical effect by means of near- 
Earth's satellites to be very careful in accounting for every tiny perturbation, of gravitational 
and non gravitational origin, which would be inevitably present in the record aliasing the 
recovered values of the quantities with which the phenomenon of interest is parameterized. 

This is the case for the Lense-Thirring drag of inertial frames [Lense and Thirring, 1918; 
Ciufolini and Wheeler, 1995] . It is a feature of the gravitational field generated by every rotating 
mass predicted by general relativity which consists in a tiny, secular precession affecting the 
node Q and the perigee uj of the orbit of a not too far freely faUing test body. For LAGEOS 
and LAGEOS 11 SLR satellites the gravitomagneic precessional rates in the field of the Earth 
amount to: 

^LAGEOS ^ simas/y, (2) 
(^LAGEOSii ^ 31.5 mas/y, (3) 
^LAGEOsn ^ -57mas/y. (4) 

The major source of uncertainty in the detection of the gravitomagnetic shift turns out to be 
the mismodeling of the first two even zonal terms of the secular classical precessions caused by 
the oblateness of the Earth. Ciufolini [1996] has proposed a strategy which could allow to detect 
the Lense-Thirring effect at 20 % precision level by using as observable a suitable combination 
of the orbital residuals of LAGEOS and LAGEOS II which would allow to cancel out all the 
static and dynamical contributions of the first two even zonal terms of the geopotential along 
with the associated errors. 

Among the various sources of systematical errors which may affect the recovery of the 
Lense-Thirring effect the Earth tides [Cartwright, 2000] play a relevant role. 

The calculations presented here, in the context of the gravitomagnetic experiment, can be 
viewed as a preliminary sensitivity test of the Lense-Thirring shift to the tidal perturbations. 
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Moreover, in view of the wide application of the two LAGEOS in many fields of space science 
they could turn out to be useful for other experiments. 

An evaluation from first principles of the amplitudes, the periods and the initial phases of 
the tidal perturbations on the nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS 
II is useful because: 

• It allows to point out which tidal constituents the Lense-Thirring shift is mainly sensitive to, 
so that people can focus the researches on them 

• It allows to get insight in how to update the orbit determination softwares. In this way the 
impact on the time series of those tidal constituents which will turn out to be more effective 
could be reduced along with the number of the parameters to be included in the least-squares 
solution 

• For a given observational time span Tabs it allows to find those constituents whose periods are 
longer than it and consequently may act like superimposed bias (e.g. the 18.6-year and 9.3-year 
zonal tides and the Ki I = 3 p = 1 q = —1 ocean diurnal tide on the perigee of LAGEOS II) 
corrupting the detection of the gravitomagnetic shift. In these cases, if we know their periods 
and initial phases we could try to fit and remove them from the time seriesQ without affecting 
also the trend of interest or, at least, it should be possible to assess the mismodeling level 
induced by them on the detection of the Lense-Thirring trend so to include also these estimates 
in the final error budget 

• It can be used as starting point for numerical simulations of the time series in order to check, 
e. g., the impact of the diurnal (m=l) and semidiurnal (m=2) tidal perturbations (not canceled 
out by the Ciufolini' s observable), as done in a further paper [Pavlis and lorio, 2001]. 

Regarding the claimed cancellation of the perturbations induced by the first two even zonal 
static and dynamical terms of the geopotential by the proposed combined residuals, it has been 
checked, as far as the tides are concerned, if it really occurs and in this case at which level of 
accuracy it takes place, and if this feature can be extended to some other tidal constituents by 
updating and extending the calculations performed with the nominal amplitudes in [Ciufolini 
et ai, 1997 pag. 2714 eq.(20), eq.(21)] 

"'^If the period of the disturbing harmonic, assumed to be known, is shorter than the time series length the 
perturbation can be rehably viewed as an empirically fit quantity. But if its period is longer than the time series 
length it is not possible to fit and remove the harmonic without removing also the true linear trend of interest. 
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The paper is organized as follows: In Sec. 2 the perturbative amplitudes on Q and uj of 
LAGEOS and LAGEOS II due to the solid tides are worked out. The calculation account 
for the frequency-dependence of Love numbers k2 and the anelasticity of the Earth's mantle. 
For some selected tidal lines -the most effective in perturbing the satellites' orbits- also the 
effect of Earth's rotation and ellipticity, accounted for by the coefficients are investigated. 
In Sec. 3 analogous calculations for the ocean tidal spectrum are performed. In Sec. 4 the 
mismodeled perturbative amplitudes of the solid and ocean tidal spectrum are compared to 
the gravitomagnetic shift over 4 years in order to check which tidal constituents, over such 
temporal span, are mismodeled at 1% level of the Lense-Thirring effect. Subsequently, the 
obtained results are employed to check extensively the feasibility of the cancellation of the 
zonal tidal constituents by means of the formula proposed by Ciufolini. Sec. 5 is devoted to 
the conclusions. 

The data employed in the present analysis are in the following table in which P is the orbital 
period and P[X] is the period of the Keplerian element X: 



dLAGEOS — 


12, 270, km 


(5) 


ttLAGEOSII = 


12, 163, km 


(6) 


^LAGEOS — 


0.0045 


(7) 


^LAGEOSII — 


0.014 


(8) 


ihAGEOS = 


110, deg 


(9) 


iLAGEOSII — 


52.65, deg 


(10) 


Plageos — 


0.1566, days 


(11) 


Plageosii = 


0.1545, days 


(12) 


P[^] LAGEOS — 


1, 043.67, days 


(13) 


P[^]lageosii — 


-569.21, days 


(14) 


P[^]lageos = 


-1,707.62, days 


(15) 


P[i^]LAGEOSII — 


821.79, days 


(16) 






(17) 
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2 Effect of the solid Earth tides on the nodal and apsidal 
lines of LAGEOS and LAGEOS II 

The free space potential [Dahlen, 1972; Smith, 1974; Wahr, 1981b; Melchior, 1983; Wang, 1997] 
by means of which the EarthQ responds to the tide generating potential $(r) [Doodson, 1921; 
Cartwright and Tayler, 1971; Cartwright and Edden, 1973; Buellsfeld, 1985; Tamura, 1987; Xi, 
1987 ; Hartmann and Wenzel, 1995; Roosbeek, 1996], at r = a for a near circular orbit satellite, 
is given by: 

= E -(— yE E Fi„U^)Gir„ie)S,mp,iuj, n,M, e). (18) 

-"^e 1=2 m=0 p=0 q=-oo 

In eq.(|l8D Fimp(i)m and Gipq{e) are the inclination and eccentricity functions respectively 
[Kaula, 1966] and the Simpq is given by: 

Si,np,{uJ, Q,M, e)= Aim E Hrkflif) cos {at + ^Jimpq - Slmf). (19) 

/ 

In eq.(|19D u, Q and are the argument of perigee, the longitude of ascending node and the 
mean anomaly, respectively, of the satellite, 9 is the Greenwich sidereal time, while Aim is given 
by: 

Aim 



2/ + l(/-m)! 
\\ 4n (Z + m)!' 

In eq.(^0]) the Condon- Shortley phase factor (—1)"* has been neglected in order to compare 
the results of this Section to those by [Bertotti and Carpino, 1989]. The phase of the generic 
harmonic in eq.(0) is constituted by three terms. The first one, at, is for every constituent 
of the tide generating potential a linear combination of the Mean Lunar Time r and the mean 
ecliptical longitudes of the astronomical motions of the Moon and the Sun [Dronkers , 1964] so 
that: 

a = 2nf = jiT + j2S + hh + j^p + j^N' + jeP^. (21) 



Since r = 6' — s eq.(21) becomes: 

(y = JiO + (j2 - Ji)s + 33h + j4P + 35N' + Ms. (22) 

^It is assumed that the Earth is endowed with a soUd inner core, a fluid outer core, and a sohd mantle 
capped by a thin continental crust, without oceans and atmosphere. Substantially, the Earth is thought as a 
set of coupled harmonic oscillators showing a variety of resonant frequencies -the normal modes- which can be 
excited by the external forcing constituents of the tide generating potential [Wahr, 1981a]. 
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The coefficients jk, k = 1, ...,Q are small integers which can assume negative, positive or null 
values. They are arranged in the so called Doodson number: 

Ji(j2 + 5)(j3 + 5).(j4 + 5)(j5 + 5)(j6 + 5), (23) 

by means of which every tidal constituent is named. In it the integer ji classifies the tides in 
long period or zonal (ji = 0), diurnal or tesseral (ji = 1) and semidiurnal or sectorial (ji = 2). 
According to the Doodson notation, the sum J2f over the tidal constituents in eg. (|19|) stands 
for the sum over the integers ji, ...je- 

The second quantity appearing in the phase of eq.(|T9|) is: 

^imp, = {I - 2p)uj + {I - 2p + q)M + m{n - 6). (24) 

The parameter: 

C^(/) = \/^f™(/)' + ^L(/f, (25) 

is the modulus of the Love number k [Love, 1926; Mathews et al., 1995] for a static, spherical 
Earth. In general, it is defined as the ratio of the free space gravitational potential 0(r), 
evaluated at the Earth' s equator, and the tide generating potential $(r) calculated at the 
mean equatorial radius. If the rotation and nonsphericity of the Earth are accounted for, there 
are different values of k for any Z,m and /. In general, they are evaluated by convolving, in 
the frequency domain, the tide generating potential with the transfer function for a non rigid 
Earth. There are several theoretical calculations for the Love numbers k [Wahr, 1981b; Wang, 
1994; Mathews et al, 1995; McCarthy, 1996; Dehant et al, 1999]. They differ in the choice of 
the Earth' s interior models [Gilbert and Dziewonski, 1975; Dziewonski and Anderson, 1981] 
adopted for the calculation of the transfer function and the departure from elasticity of the 
mantle' s behavior. Particularly important is also if they account for, or not, the normal modes 
in the diurnal band. In eq.(^) the lERS conventions [McCarthy, 1996] have been used for 

kUf) = Imh^ + 5kf^\ (26) 
kll{f) = Reki„, + 6kf, (27) 
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where Im kim and Re ki^ are the frequency- independent parts of Love numbers, and 6k j , 
6k'j are the frequency-dependent corrections. They are important in the diurnal band, through 
(J/cj , and in the zonal band with dk'p^^ due to the anelasticity of the mantle. 

The factor 5imf is the phase lag of the response of the solid Earth with respect to the tidal 
constituent of degree /, order m and circular frequency a induced by the anelasticity in the 
mantle [Varga, 1998]: 

tan5z™/ = ^. (28) 

Note that if kj^^f) = 0, also tanSimf = 0. It may be important to know these parameters. 
Indeed, the initial phases of some semisecular tidal constituents, like the 18.6-year tide, consist 
entirely of them and the astronomical longitudes of the Sun and the Moon. Indeed eq.(P^) 
vanishes for the even {I — 2p = 0) zonal (m = 0) long period {I — 2p + q = 0) perturbations. 
For a time span Tots shorter than their periods they could be fitted and removed from the time 
series provided that their periods and initial phases are exactly known. The same remarks will 
hold for the ocean tides. 

The quantities Hf^{f) appearing in eq.(0) are the coefficients of the harmonic expansion of 
the tide generating potential. They contain the lunisolar ephemerides information and define 
the modulus of the vertical shift V/g in the equipotential surface at the Earth's surface for 
r = with respect to the case in which only proper Earth's gravity is considered; g is the 
acceleration of gravity at the equator. The values used in the present calculation are those 
recently released by Roosbeek [1996]. Their accuracy is of the order of 10~^ m. In order to 
use them in place of the coefficients of Cartwright and Edden [1973], which have a different 
normalization, a multiplication for a suitable conversion factor [McCarthy, 1996] has been 
performed. 

The eg. ([T8|) is the dynamical, non central part of the geopotential due to the response of 
the solid Earth to the forcing lunisolar tidal solicitation. In the linear Lagrange equations of 
perturbation theory for the rates of change of the Keplerian elements of a satellite it plays the 
role of the perturbative function TZ. The latter is the non central part of the total mechanical 
energy of the satellite, with the sign reversed, according to the geophysical convention. The 
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Lagrange equations are [Kaula, 1966]: 

da 2 dn 



dt na dM' 

de l-e2 on fl- 6^)1/2 



dt na^e dA4 na?e du ' 

di _ . 1 dH _ 

dt na^{l — e^y/'^sm.i du 

1 &R 



na^(l — e^)^/^ sini dVL^ 

dn _ 1 dn 

dt na^(l — e^)^/^ sini di 

du _ . 1 &R, 

dt na^ {1 — e^y/^ sin i di 



(29) 
(30) 

(31) 
(32) 



(1 - f 2)1/2 on , , 

+^ ^ IT' 33 

dM _ l-e^ dn ^ dn 
dt na^e de na da 

In it n = {GM^y/'^a''^^'^ — 2t^/P is the mean motion of the sateUite. It is an observed fact that 
the secular motions are the dominant perturbation in the elements a;, of geodetically 

useful satellites. So, taking as constants a, e, i and consider as linearly variable in time a;, Q, M. 
and 6*, besides at, the following expressions, at first order, may be worked out: 



n oo / n 



na 



^Vl — sini 



X 



=0 m=0 



xA^E E -^GivA^tlHr^^n^nrnv,, (35) 



p=0 g=— oo •'P 



O oo « n 



na- 



=71 



X 



=0 m=0 



where: 



v4 F ^Gipi cosi dFimp ^ 1 ^ i-W ff™ dr, ^ ^-^fi^ 

p=0 q=-oo "'^ ' 



= (/ - 2p)a; + (/ - 2p + g)M + m(0 - 9) + at - Simf, (37) 
fp^{l- 2p)Cj + {l-2p + q)M + m(Q - ^) + (7. (38) 
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Only the perturbations whose periods are much longer than those of the orbital satellites 
motions, which, typically, amount to a few hours, are relevant for many geophysical and astro- 
nomical applications. This implies that in the summations of eqs. (^5D -(|36|) only those terms in 
which the rate of the mean anomaly does not appear must be retained, i. e. the condition: 

l-2p + q = (39) 

must be fulfilled. Moreover, if the effect of Earth's diurnal rotation, which could introduce 
periodicities of the order of 24 hours, is to be neglected, one must retain only those terms 
in which the non negative multiplicative coefficient ji of the Greenwich sidereal time in at 
coincides to the order m of the tidal constituent considered: in this way in fp the contributions 
of 9 are equal and opposite, and cancel out. With these bounds on I, m, p and q the circular 
frequencies of the perturbations of interest become: 

= + (/ _ 2p)i:u + rati (40) 

in which: 

T/ = (j2 - rn)s + js/i + j^p + j^N' + jePs- (41) 

In the performed calculation only the degree / = 2 constituents have been considered due to 
the smallness of the and H^{f). For I = 2, p runs from to 2, and so, in virtue of the 
condition I ~ 2p + q = 0, q assumes the values —2, 0, 2. From an inspection of the table of 
the eccentricity functions Gipq{e) in [Kaula, 1966] it turns out that G20-2 = ^222 = 0, while 
^210 = (1 ~ e^)^^/^. For this combination of /, p and q the relation / — 2p = is fulfilled: the 
frequencies of the perturbations are, in this case, given by: 

fp = tf + mh. (42) 

In Tab.[^, Tab.^ and Tab.^ the results for the nodal lines of LAGEOS and LAGEOS II, 
and the apsidal line of LAGEOS II are shown; since the observable quantity for u is eau 
Ciufolini, 1996|| , the calculation for the perigee of LAGEOS, due to the notable smallness of 



the eccentricity of its orbit, have not been performed. The tidal lines for which the analysis 
was performed have been chosen also in order to make a comparison with the perturbative 
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amplitudes of the ocean tides based on the results of EGM96 gravity model which will be 
shown in the next Section. 

The results presented in Tab.|l], Tab.| and Tab.|| can be compared to those of Dow [1988] 
and Carpino [Bertotti and Carpino, 1989]. In doing so it must be kept in mind that both 
these authors have neglected not only the contribution of k'^ but also the anelasticity and the 
frequency dependence of the Love numbers k^^. Moreover, Dow includes in his analysis, in a not 
entirely clear manner, also the indirect influences of the oblateness of the Earth [Balmino, 1974; 
Dow, 1988; Casotto, 1989]. Obviously, in their analyses the LAGEOS II is not present since it 
was launched only in 1992. Another important factor to be considered is the actual sensitivity 
in measurements of Q and u, in the sense that the eventual discrepancies between the present 
results and the other ones must be not smaller than the experimental error in the Keplerian 
elements if one wishes to check the theoretical assumptions behind the different models adopted. 
Carpino has analyzed the inclination and the node of LAGEOS only. His value for the zonal 
18.6-year tide is —1087.24 mas, while Tab.|I| gives —1079.38 mas; the difference amounts to 7.86 
mas, the 0.72 % of the "elastic", frequency-independent Carpino's value. Considering that for 
the node Q of LAGEOS the present accuracy is of the order of the mas, 7.86 mas could be in 
principle detected, allowing for a discrimination between the different models adopted in the 
calculation. For the Ki tide, one of the most powerful constraints in perturbing the satellites' 
orbits, Tab.|l] quotes 1744.38 mas against 2144.46 mas of Carpino's result; the gap is 400.08 mas, 
the 18.6 % of Carpino. In the sectorial band, the present analysis quotes for the K2 —92.37 
mas while Carpino has —97.54 mas; there is a difference of 5.17 mas, the 5.3 % of the Carpino's 
value. It must pointed out that there are other tidal lines for which the difference falls below 
the mas level, as is the case for the 9.3-year tide. As it could be expected, the major differences 
between the present "anelastic" , frequency-dependent calculations and the other ones based on 
a single, real value for the Love number k2 lie in the diurnal band: in it the contribution of 
anelasticity is not particularly relevant, but, as already pointed out, the elastic part of k2i is 
strongly dependent on frequencies of the tidal spectrum. May be interesting to note that when 
the calculation have been repeated with the same value ^2 = 0.317 adopted by Carpino, his 
results have been obtained again. 

Up to now the effects induced by the Earth's flattening and the Earth's rotation have been 
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neglected. If they have to be analyzed, it is necessary to take in the free space potential' s 
expansion the part [Wahr, 1981b]: 

oo I p 

Regi: T.T.HrU){-r'kUf) YUM. A) (43) 

1=2 m=0 f 

and work out it in the same manner as done for the spherical nonrotating Earth's contribution. 
In applying the standard transformation to the orbital elements of [Kaula, 1966] one has to 
substitute everywhere I with 1 + 2. So the equations for the perturbations on the node and the 
perigee become: 

n oo i p 

na"^V 1 — sm« r 

I +00 

X^imyE J2 "'^ g/+2 pq Sin7fl+2 mpq, (44) 

Jp p=Q q=-oo 

oo / 



oo J p -| 

1=0 m=0 ' fp 
1+2 pq 



1+2 mp' 



p=0 q=-oo ^ 

COSidFi+2 mp^ 1 . f.^. 

~ ^ di sm7j,+2 ^pg. (45) 

The corrections induced by the Earth's flattening and rotation, due to the smallness of k^, 
have been calculated only for those tidal lines which have resulted to be the most powerful 
in perturbing the node and the perigee of LAGEOS and LAGEOS II, i.e. the zonal 18.6-year 
tide and the Ki. Adopting the values quoted in [Dehant et al, 1999] for k'^ and the eq.(^) 
and eq.(^5[) it is possible to obtain the results summarized in Tab.^. The quoted values, with 
the exception of the node of LAGEOS II, fall below the mas level resulting, at the present, 
undetectable. 

Concerning the mismodeling of the solid tidal perturbations of the nodes of LAGEOS and 
LAGEOS II and the perigee of LAGEOS II, particularly important for the gravitomagnetic 
experiment, by inspecting the analytical expressions of their perturbations given in eqs.(^)- 
( PB| ) it is possible to note that almost all quantities entering in them are very well determined 
with the possible exception of the inclination i and the Love number k. Concerning the in- 



clination i, by assuming 6i = 0.5 mas [Ciufolini, 1989], we have calculated 



dA{Q) 
di 



6i and 
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Si for the 18.6-year tide, Ki, and 5*2 which are the most powerful tidal constituents 
in perturbing LAGEOS and LAGEOS II orbits, as can be inferred from Tab.^Tab.|^. The 
results are of the order of 10~^ mas, so that we can neglect the effect of uncertainties in the 
inclination determination. About the Love numbers k2, we have assessed the uncertainties in 
them by calculating for certain tidal constituents the factor 5/22/^2; ^2 is the average on the 
values released by the most reliable models and Sk2 is its standard deviation. According to the 
recommendations of the Working Group of Theoretical Tidal Model of the Earth Tide Com- 
mission ( littp: / /www. astro. oma.be/Dl/EARTH_TIDES/wgtide. htmlp , in the diurnal band we 
have chosen the values released by Mathews et al. [1995], McCarthy [1996] and the two sets 
by Dehant et al. [1999]. For the zonal and sectorial bands we have included also the results of 
Wang [1994]. The uncertainties calculated in the Love numbers k2 span from 0.4% to 1.5% for 
the tides of interest. 

3 Ocean tidal perturbations on the nodal and apsidal 
lines of LAGEOS and LAGEOS II 

The effects of the ocean tides [Defant, 1961; Neumann et al, 1966; Pekeris and Akkad, 1969; 
Hendershott and Munk, 1970; Hendershott, 1972; 1973; Schwiderski, 1980; 1981; Gill, 1982, 
Zahel, 1997], for a constituent of given frequency /, can globally accounted for by means of the 
potential: 



1=0 m=0 + ' p=0(?=-oo 



COS 7^;, 



± T l—m even 
mpq 



i 

mpq 



(46) 



l—m odd 



in which: 



^timpq = - 2p)^ + {l-2p + q)M + m{n - 6) ± {at - etj), (47) 



In eq.(^6]) expressions like X)+ cos (a ± b) mean A'^ cos {a + b) + A" cos {a — b). The sign 
-|- refers to the progressive (westwards) waves while the sign — is for the waves moving from 
West to East. Moreover, p^, denotes the water density [McCarthy, 1996], are the load Love 
numbers [Farrell, 1972] accounting for the ocean loading [Jentzsch, 1997], the Cj^j are the 
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ocean tidal coefficients of EGM96 gravity model [Lemoine et. al, 1998] and are the phase 
shifts due to hydrodynamics of the oceans. 

The equations for the tidal ocean perturbations may be worked out as already done for the 
solid Earth tides in the previous Section. At ffist order, one obtains: 

1 oo / 



R 



•X 



na^y I — Sim i^Q^^Q _,_ r 

-1 l—m even 

(49) 



I +00 j-p -1 



Imp 



sill Iflmpq 



p=0 q=—oo 

oo I 



COS 7^;, 



± 

mpq 



l—m odd 



1 oo t — p 

na V i e + ' 



X 



I +00 1 _ „2 J/^ nncorlFP 1 r oivi -fl-m even 

E\ ^ r -L e Uj^ipq L.UD t tirif^ip , 1 



p=0 g=— oo 



-COS 7^;, 



± 

mpq 



(50) 

/— m odd 



in which: 



fp = {l- 2p)dj + {l-2p + q)M + m{n -e)±cx. (51) 

It should be noted that the frequencies of the perturbations given by the eq.(^) are different, 
in general, from the frequencies of the solid Earth tidal perturbations given by eq.(^). While 
for the solid tides the diurnal modulation due to 6 cancels out automatically if one considers 
those terms in which ji = m, for the ocean tides, in general, this does not happen because 
of the presence of the Eastwards waves due to the non equilibrium pattern of the ocean tidal 
bulge. Long periodicities can be obtained considering those combinations of I, p and q for 
which l — 2p + q = and retaining the Westward prograde terms with ji = m. Only in this way 
in eq.(|5l|) the contributions of 9 are equal and opposite, and cancel out. With these bounds on 
/, m, p and q the frequencies of the perturbations of interest become: 

/p = f / + (/ - 2p)(:o + mh. (52) 

It is worthwhile noting that the frequencies fp are identical to those of solid tidal perturbations, 
so that the satellites cannot distinguish one effect from the other. For / = 2, as for the solid 
tides, I — 2p = holds. The frequencies of the perturbations are, in this case, given by: 

fp = tf + mtl. (53) 
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If / is odd the situation is different because now all the terms with p running from to I must 
be considered and q is different from zero: q = 2p — l. So, the rates of the perturbations include 
the contribute of u. 

The eqs.(^9|)-(|50|) have been adopted in order to compute the amplitudes of the ocean tidal 
perturbations on the Keplerian elements Q, u, M. and i for LAGEOS and LAGEOS II. The 
calculation have been performed, considering only the progressive waves, for the following tidal 
lines: 

Mm (065.455), Sa (056.554), Mf (075.555), S^a (057.555); 

Ki (165.555), Oi (145.555), Pi (163.555), Qi (135.655); 

K2 (275.555), M2 (255.555), S2 (273.555), N2 (245.655), T2 (272.556). 

For each of these tidal constituents the following terms have been calculated: 

/ = 2, p = 1, g = because the eccentricity functions Gipg{e) for p = 0, q = —2 and 

p = 2, q = 2 vanish. 

/ = 3, p = 1, q = —1 and 1 = 3, p = 2, q=l because G30-3 and G333 are not quoted in (Kaula 
1966) due to their smallness: indeed, the Gipq{e) are proportional to e''^'. 

I = 4, p = 2, q = because the other admissible combinations of /, p and q give rise to negligible 
eccentricity functions. So, also for / = 4 the condition / — 2p = 0, in practice, holds and the 
constituents of degree / = 2 and I = 4 generate detectable perturbations with identical periods. 
Similar analyses can be found in [Christodoulidis, 1978; Dow, 1988]. When Christodoulidis 
performed his study, which is relative to only five constituents, neither the LAGEOS nor the 
LAGEOS II were in orbit, while Dow has sampled the tidal spectrum for LAGEOS in a poorer 
manner with respect to this study in the sense that, for each constituent, only the terms of 
degree / = 2 have been considered with the exception of the Ki whose / = 4 contribution has 
been also analyzed. Moreover, when these work have been realized there were a few coefficients 
Qm/ available with relevant associated errors. 

It should be pointed out that the values obtained for the coefficients Cj^j are, in general, 
biased by the effects of the anelasticity of the solid Earth's mantle and by all other phenomena 
which have not been explicitly modeled in the nominal background models used [Lemoine et 
ai, 1998]. For example, in the values recovered for the Sa are included climatological effects 
which have not gravitational origin; in the S2 coefficients are also included the variations of the 
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atmospheric pressure due to the atmospheric tides. 

In Tab.^ Tab.|^ and Tab.^ the present resuhs for the nodes Q of LAGEOS and LAGEOS 
II, and the argument of perigee u for LAGEOS II are quoted. 

From an accurate inspection of TabJ^, Tab.0 it is possible to note that, for the nodes fl, 
only the even degree terms give an appreciable contribute; the / = 3 part of the spectrum 
is totally negligible. This so because Afij is proportional to the Gipq{e) functions which, 
in turn, are, in general, proportional to e'"^': in this case the eccentricity functions used are 

G(e)3i-i=e(l-e2)-5/2 = G(e)32i. 

Among the long period zonal tides, the Solar annual tide Sa (056.554) exerts the most 
relevant action on the nodes, with an associated percent error in the amplitudes of 6.7 %. 
It is interesting to compare for this tidal line the ocean tidal perturbations of degree / = 2 
^LAGEOsi^) = -20.55 mas, AfX'cEosiii^) = 37.7 mas, with those due to the sohd Earth tides 
^LAGEOsi^) = 9-96 mas, Afj^^Eosiii^) = -18.28 mas. The ocean amplitudes amount to 204 
% of the solid ones, while for the other zonal constituents they vary from 9.9 % for the to 
the 18.5 % of the Ssa- This seems to point toward that the recovered value of 6*2^^ for the Sa is 
biased by other climatological effects than the tides; indeed, for all the other tidal lines, zonal or 
not, the ocean tidal perturbations of degree / = 2 amount to almost 10 % of the corresponding 
solid Earth tides perturbations. 

The terms of degree / = 2 of the tesseral tides Ki (165.555) and Pi (163.555) induce very 
large perturbations on the node of LAGEOS and, to a lesser extent, of LAGEOS II; TabJ^ and 
Tab.|^ quote 156.55 mas, —11.49 mas for the former and —35.69 mas, —3 mas for the latter. 
The associated percent errors are 3.8 % and 8.1 %, respectively, for the two constituents. For 
the terms of degree I = 4, which have the same periods of those of degree / = 2, the situation 
is analogous: Tab.|^ and Tab.|^ quote 4.63 mas and —0.32 mas for LAGEOS and 41.58 and 0.7 
mas for LAGEOS II. The percent errors associated with Ki and Pi, for / = 4, are 3.9 % and 
8 % respectively. For all the tesseral lines investigated the ocean tidal perturbations of degree 
/ = 2 are, in general, the 10 % of the solid tidal perturbations. 

Among the sectorial tides, the most relevant in perturbing the nodes of LAGEOS satellites, 
on large temporal scales, is the 5*2 (273.555): Tab.^and Tab.|^ quote 9.45 mas (/ = 2) and 15.08 
mas (/ = 4) for LAGEOS, and 6.87 mas (/ = 2) and 6.79 mas (/ = 4) for LAGEOS II. The 
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associated percent errors are 3.9 % for / = 2 and 5.2 % for I = 4. The 5*2 ocean perturbation 
of degree / = 2 amounts to nearly 5 % of the corresponding sohd tide. 

About the zonal 18.6-year and 9.3-year tides, recently both Starlette and LAGEOS SLR 
satellites passed their 19th year in orbit and this span of time is now adequate to get reliable 
information about these tides [Cheng, 1995; Eanes, 1995] due to their slow frequencies they can 
be adequately modeled in terms of the equilibrium theory through the Hj^ coefficients and a 
complex Love number accounting for the anelasticity of the mantle. So, concerning them the 
results quoted for the solid tides can be considered adequately representative. 

In Tab.^] the amplitudes of the perturbations on the argument of perigee u for LAGEOS II 
are quoted. For this orbital element the factor: 

l-e^ dGipg cosidFimp^ 

^imp — ; : — : — r — '-^ipq [^^) 

e ae smi ai 

to which Aujf is proportional makes the contributions of the / = 3 terms not negligible. For the 
even degree terms the situation is quite similar to that of Q in the sense that the most influent 
tidal lines are the Sa, Ki, Pi and 5*2. 

Once again, among the long period tides the Sa exhibits a characteristic behavior. Indeed, its 
/ = 3 contributions are much stronger than those of the other zonal tides. This fact could be 
connected to the large values obtained in its C^qj- coefficient and people believe that it partially 
represents north to south hemisphere mass transport effects with an annual cycle nontidal in 
origin. The / = 3 terms present, in general, for the perigee of LAGEOS II a very interesting 
spectrum also for the tesseral and sectorial bands: for LAGEOS II there are lots of tidal lines 
which induce, on large temporal scales, very relevant perturbations on u, with periods of the 
order of an year or more. Above all, it must be quoted the effect of Ki line for p = 1, q = — 1: 
the perturbation induced amounts to -1136 mas with period of -1851.9 days. These values are 
comparable to the effects induced by the solid Earth tidal constituent of degree Z = 2 on the 
node Q of LAGEOS. By comparing the / = 2 terms with the corresponding solid tides, it can 
be noted that also for the perigee the proportions are the same already seen for the nodes. 

About the mismodeling of the ocean tidal orbital perturbations, which are the worst de- 
termined part of the spectrum of the Earth tidal response, the major source of uncertainties 
in their amplitudes resides in the EGM96 coefficients Cjj^f and the ocean loading parameters 
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k'l. Concerning the ocean loading, Pagiatakis [1990] in a first step has recalculated k'l for an 
elastic, isotropic and non-rotating Earth: for / < 800 he claims that his estimates differ from 
those by Farrell [1972], calculated with the same hypotheses, of less than 1%. Subsequently, 
he added to the equations, one at a time, the effects of anisotropy, rotation and dissipation; 
for low values of / their effects on the results of the calculations amount to less than 1%. By 
inspecting Tab.B|-Tab.|] it has been decided to calculate ^fc; of the perigee of LAGEOS 

II for Ki I = "i p = 1. First, we have calculated mean and standard deviation of the values 
for fcg by Farrell and Pagiatakis obtaining 5k'^/k'^ = 0.9%, in according to the estimates by 
Pagiatakis. Then, by assuming in a pessimistic way that the global effect of the departures 
from these symmetric models yield to a total 5k'^/k'-^ = 2%, we have obtained 5uj^^ = 5.5 mas. 
Subsequently, for this constituent and for all other tidal lines we have calculated the effect of 
the mismodeling of Cj^j- as quoted in EGM96. 

4 The influence of the tides on the detection of the 
Lense-Thirring drag 

The theoretical general relativistic expressions for the gravitomagnetic precessions of the node 
and the perigee of a test particle in the field of a massive rotating central body -the Lense- 
Thirring effect- are given by: 

= ia3(l-^e2)3/2' 
G 2J 

^LT = -^ ^3(i_e2)3/2 3cQsz, (56) 

in which: 

J angular momentum of the central rotating body, g CTn?s~^. For the Earth its value is 

5.9 ■ 10^° g cm'^s^\ 

c speed of light in vacuo, cm s^^. 

According to [Ciufolini, 1996; Ciufolini et al. 1997; 1998], it should be possible to detect 
the Lense-Thirring shift at a 20 % level of accuracy through the formula: 

5tl^ + ci5tl^^ + C256j^^ = 60.05 x ixlt, (57) 
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where Ci ~ 0.295, C2 — —0.35, fi^x is an adimensional scale parameter which is 1 in general 
relativity and in Newtonian mechanics and , 6Q^^ , 6uj^^ are the residuals, in mas, of the 
nodes of LAGEOS and LAGEOS 11 and the perigee of LAGEOS 11 calculated, e.g., with the 
aid of the NASA-Goddard software GEODYN 11. In the analyses performed with the real data 
[Ciufolini et al, 1997; 1998] the Lense-Thirring parameter /i has been purposely left out the 
dynamical models of GEODYN 11 by setting = so that the residuals can entirely account 
for it. General relativity predicts for eq.(^) a secular trend with a slope of 60.05 mas/y if it 
is calculated by means of eqs.(^5])-(0). Eg. (|57|) has been recently investigated in [Ciufolini et 
ai, 1998] for a 4 years time span. 

In view of a refinement of the error budget of this experiment the results obtained in Sec. 
2 and Sec. 3 have been used in order to calculate the mismodeled amplitudes of the solid and 
ocean tidal perturbations , and 6uj^^; they have been subsequently compared with the 
gravitomagnetic precessions over 4 years Afi^^ = 124 mas, Afi^^ = 126 mas and Auoi^T = -228 
mas. In Tab.^we have quoted those solid tidal lines whose mismodeled perturbative amplitudes 
amount to 1%, at least, of the gravitomagnetic shifts. It turns out that only the 18.6-year tide 
and Ki exceed this cutoff. Regarding the ocean tides, it turns out from Tab.^ that the perigee 
of LAGEOS 11 is more sensitive to the mismodeling of the ocean part of the Earth response 
to the tide generating potential. In particular, the effect of Ki I = 3 p = 1 q = —1 is relevant 



dAjuj) 

act , 



dA{u,) 
dk' 



6k[ of 64.5 mas amounting to 28.3 % of Ault over 4 



with a total 5u 
years. 

The combination i/lt is useful since it should vanish if calculated for the even zonal con- 
tributions C20 and C40 of the geopotential [Ciufolini, 1996]. More precisely, the right side of 
eq.(|57D should become equal to zero if the left side were calculated for any of these two even 
zonal contributions, both of static and dynamical origin; the nearer to zero is the right side, 
the smaller is the systematic uncertainty in fi due to the contribution considered. 

In order to test preliminarily this important feature for the case of tides, in a very conser- 
vative way the results obtained in the present work have been used in eq.(0) assuming, for 
the sake of clarity and in order to make easier the comparison with [Ciufolini et ai, 1997], a 
1 year time span and the nominal values of the calculated tidal perturbative amplitudes, as if 
the zonal solid and ocean tides were not at all included in the GEODYN 11 djTiamical models 
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so that the residuals should account entirely for them. The results are released in Tab.p!0| and 
Tab.|TT]. The values of 5/i quoted there for the various zonal tidal lines may be considered as the 
systematic error in fi due to the chosen constituents, if considered one at a time by neglecting 
any possible reciprocal correlation among the other tidal lines. Tab.|lO| and Tab.pT] show that 
the percent error in the general relativistic value of /i due to the 18.6-year tide, the most dan- 
gerous one in recovering the LT since it superimposes to the gravitomagnetic trend over time 
spans of a few year, amounts to 21.9 %, while for all the other zonal tides it reduces to 0.1 % 
or less. This means that, even if fully present in the combined residuals, the Z = 2 m = tides, 
with the exception of the 18.6-year tide, do not affect the recovery of by means of i/lt so 
that there is no need including them into the final least-squares scheme. It is interesting to 
compare the present results to those released in [Ciufolini et al. 1997] for the 18.6-year tide. 
The value -0.219 due to the solid component for /i quoted in Tab.|TO| must be compared to 
-0.361 in [Ciufolini et al. 1997], with an improvement of 39.3 %. In the cited work there is no 
reference to any estimate of the mismodeling of the 18.6-year tide, so that we have used the 
nominal tidal perturbative amplitudes released in it: AiVt^) = —997 mas, A{Q^^) = 1805 mas 
and A{uj^^) = —1265 mas. These figures for the perturbative amplitudes due to the solid Earth 
tide of 18.6-year are notably different from those quoted in the present study. In [Ciufolini et 
al., 1997] the theoretical framework in which those numbers have been calculated (F. Vespe, 
private communication, 1999) is based on the assumption of a spherical, static, elastic Earth 
with a single nominal value of k2 = 0.317 used for the entire tidal spectrum. The inclusion of 
the tiny corrections due to the Earth's flattening and rotation on the perturbative amplitudes 
of Q and u for the 18.6-year tide could allow to slightly improve the related uncertainty in fi; 
it would amount to 20.6 %. But since the present-day accuracy in laser ranging measurements 
could hardly allow to detect these small effects, their utilization in eq.(^) is debatable. 

Remember that the result quoted for the 18.6-year tide is obtained in the worst possible 
case, i.e. a time span of only 1 year and the assumption that the residuals have been built up 
by neglecting completely the zonal tides in the dynamical models. Recall that if the residuals 
accounted for the 18.6-year tide it would not be possible to view it as an empirically fit quantity 
unless a Tabs of, at least, 18.6 years is adopted. However, if, more realistically, we calculate 
eq. (|57D with the mismodeled amplitudes quoted in the first row of Tab.^for the 18.6-year tide 
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we obtain 5^ =< —3.51 • 10^'^. This conclusion has been confirmed also by calculating the time 
average of the 18.6-year tide over different time spans by using the mismodeling level adopted 
here [Pavlis and lorio, 2001]. 

Even though a cancellation is not expected as for the first two even zonal constituents, 
calculating the left hand side of eq.(|^) for the other tides yields, at least, an order of magnitude 
of their effect on /i. An interesting, unpredicted feature stands out for the odd zonal ocean 
tides. The contribution of / = 3 zonal ocean tidal nominal perturbations over 1 year to 6^ can 



be found in Tab.^ and Tab.|l^. From an inspection of them it is clear that the sensitivity of 
perigee of LAGEOS II to the / = 3 part of the ocean tidal spectrum may affect the recovery 
of the Lense-Thirring parameter /i by means of the proposed combined residuals, especially 
as far as Sa and Ssa are concerned. This fact agrees with the results of Tab.^ which tells us 
that the mismodeled parts of Sa and Sga are not negligible fractions of Auj[^. However, if the 
mismodeled amplitudes are employed in eq.(p7D it can be seen that, over 1 year, a cancellation 
of the order of 10^^ {Sa I = 3 p = 2) and 10^^ {Sa I = 3 p = 1; Sga I = 3 p = 1,2) takes place. 
The contributions of the mismodeling on Mm and Mf are completely negligible. So, we can 
conclude that also the / = 3 part of the zonal ocean tidal spectrum may not affect the combined 
residuals in a sensible manner if the / = 3 part of Sa and Ssa is properly accounted for. 

5 Conclusions 

The detection of the general relativistic Lense-Thirring drag in the field of the Earth by means 
of the combined residuals of the two LAGEOS laser-ranged satellites is affected, among other 
factors, also by the Earth solid and ocean tidal perturbations. 

Concerning the solid tides of degree 1 = 2, the most effective constituents turn out to be 
the semisecular 18.6-year, the Ki and the 5*2 which induce on the examined orbital elements 
perturbations of the order of 10^ — 10'^ mas with periods ranging from 111.24 days for the 5*2 
on LAGEOS II to 6798.38 days for the 18.6-year tide. 

Regarding the ocean tides, we have analyzed the degree / = 2, 3, 4 terms of 13 constituents 
by using the data of the EGM 96 model. If, from one hand, the nominal ocean amplitudes of 
degree I = 2 amount to almost 10% of the solid ones for the same degree, from the other hand 
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the obtained results for the perigee of LAGEOS II show that this orbital element is very sensitive 
to the degree / = 3 part of the tidal spectrum. Indeed, for many tidal lines the perturbations on 
it with I = 3 p = 1,2 q = —1,1 exhibit long periods of some years and amplitudes of the order 
of magnitude of those due to the sohd Earth tides for the node. The Kil — 3p — Iq — —1 
line stands out with its period of -1851.9 days (5.07 years) and amplitude of 1136 mas. Such 
term has a 5.2 % error. 

The calculations performed here have been used, in view of a refinement of the error budget 
of the gravitomagnetic LAGEOS experiment, in order to check preliminarily which tidal con- 
stituents arc really important in perturbing the combined residuals so to fit and remove them 
from the data, if possible, or, at least, to evaluate the systematic error induced by them. Tab. 
5 and Tab. 9 show that, over a 4 years time span the nodes of the two LAGEOS are sensitive to 
the even components of the 18.6-year fine, the Ki and the Sa at a 1% level at least. Moreover, 
the perigee of LAGEOS II turns out to be very sensitive to the I — 3 part of the ocean tidal 
spectrum. 

Concerning the very long period tides, an accurate calculation from first principles of their 
amplitudes, periods and initial phases turns out to be important in the case of an observational 
time span shorter than their periods. Indeed, while for a biasing true linear trend one can only 
assess as more accurately as possible the error induced by it, in the case of the arc of sinusoid 
of a harmonic perturbation whose period is longer than the adopted time span it should be 
possible to fit and remove it from the signal, without affecting the trend of interest, provided 
that its period and initial phases are exactly known. In any case it could be possible to assess 
reliably its role in the error budget. 

However, in the Lense-Thirring experiment we have shown in a preliminary way that this 
should not be necessary, at least for the semisecular even zonal 18.6-year tide independently 
of the time series length. Indeed, if calculated at the level of accuracy shown in this paper, 
it affects the combined residuals at a level < 10~^. The other even zonal tides do not create 
problems. Also the / = 3 m = tides, and this is an unpredicted feature cancel out at a 
level < 10~^-10~^ if they are properly accounted for in building up the residuals. The results 
presented in this paper not only confirm the usefulness of the formula by Ciufolini in canceling 
out the I — 2,4: m — tides, but also extend its validity to the I — 3 m — part of the tidal 
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response spectrum. 

The impact of the tesseral (m = 1) and sectorial (m = 2) tidal constituents on the combined 
residuals, which are fully sensitive to them, over different time spans and time steps will be the 
subject of a forthcoming numerical analysis based on the results presented here. 
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Solid Earth tidal perturbations on the node of LAGEOS 
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Table 1: Pcrturbativc amplitudes on the node fl of LAGEOS due to solid Earth tides. In the first column the 
Doodson number of each constituent is quoted followed by the Darwin's name, when it is present. The tidal 
lines are listed in order of decreasing periods. The coefScients Hp{f) are those recently calculated by Roosbeek 
and multiplied by suitable normalization factors (lERS standards) in order to make possible a comparison with 
those of Cartwright and Edden. tanSimf expresses the phase lag of the solid Earth response with respect to 
the tidal potential due to the anelasticity in the mantle. 
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Solid Earth tidal perturbations on the node fl of LAGEOS II 
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Table 2: Perturbative amplitudes on the node f2 of LAGEOS II due to solid Earth tides. In the first column 
the Doodson number of each constituent is quoted followed by the Darwin's name, when it is present. The tidal 
lines are listed in order of decreasing periods. The coefficients Hp{ f) are those recently calculated by Roosbcek 
and multiplied by suitable normalization factors (lERS standards) in order to make possible a comparison with 
those of Cartwright and Edden. ta.nSimf expresses the phase lag of the solid Earth response with respect to 
the tidal potential due to the anelasticity in the mantle. 



31 



Solid Earth tidal perturbations on the perigee uj of LAGEOS II 


Tide 


1=2, p=l, q=0 


Period (days) 


Amplitude (mas) 


^2 Love number 


(meters) 


tifi.Tl rt/™ f 
ijUiJ-x \JlffiJ 


055 565 

yj \j \j * yj ^-j 


6798 38 


-1375 58 


0.315 


02792 


-0.01715 


055.575 


3399.19 


6. 66 


0.313 


-0.000272 


-0.015584 


05().-')-')4 S'„ 


.3(i5.27 


12. (i9 


0.307 


-0.00492 


-0.01135 


057.555 Ssa 


182.62 


39.77 


0.305 


-0.03099 


-0.01029 


065.455 M„ 


27.55 


6.74 


0.302 


-0.03518 


-0.00782 


075.555 Mf 


13.66 


6.30 


0.301 


-0.06659 


-0.007059 


165.565 


-621.22 


290.43 


0.254 


0.050028 


-0.005653 


165.555 Ki 


-569.21 


1982.14 


0.257 


0.3687012 


-0.0055933 


165 545 


-525.23 


-36.52 


0.259 


-0.007295 


-0.005541 


163.555 Pi 


-138.26 


-177.56 


0.286 


-0.1219 


-0.005017 


145.555 Oi 


-13.33 


-38.16 


0.297 


-0.26214 


-0.00484 


135.655 Qi 


-8.98 


-4.92 


0.297 


-0.05019 


-0.00483 


275.555 K2 


-284.6 


-88.19 


0.301 


0.079915 


-0.004318 


274.556 


-159.96 


-0.38 


0.301 


0.000625 


-0.00431 


274.554 


-159.95 


1.52 


0.301 


-0.00246 


-0.0043 


273.555 52 


-111.24 


-126.83 


0.301 


0.29402 


-0.00431 


272.556 T2 


-85.27 


-5.68 


0.301 


0.017188 


-0.00431 


255.555 M2 


-13.03 


-31.9 


0.301 


0.6319 


-0.004318 


245.655 N2 


-8.84 


-4.15 


0.301 


0.12099 


-0.00431 



Table 3: Perturbative amplitudes on the perigee w of LAGEOS II due to solid Earth tides. In the first column 
the Doodson number of each constituent is quoted followed by the Darwin's name, when it is present. The tidal 
lines arc listed in order of decreasing periods. The coefficients HJ^{f) are those recently calculated by Roosbeek 
and multiplied by suitable normalization factors (lERS standards) in order to make possible a comparison with 
those of Cartwright and Edden. tan Simf expresses the phase lag of the solid Earth response with respect to 
the tidal potential due to the anelasticity in the mantle. 



Corrections due to Earth's flattening and rotation to solid tidal perturbations 
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Table 4: Corrections to the Earth solid tidal perturbations on the nodes of LAGEOS and the perigee of 
LAGEOS II due to the Earth's flattening and Earth's rotation. The values adopted for are those quoted 
by Dehant [1999]. 



Mismodeled solid tidal perturbations on nodes of LAGEOS and LAGEOS II and perigee of LAGEOS II 
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Table 5: Mismodeled solid tidal perturbations on nodes il of LAGEOS and LAGEOS II and the perigee u) of 
LAGEOS II compared to their gravitomagnetic precessions over 4 years. 
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Ocean tidal perturbations on the node S2 of LACiiEOS 
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Table 6: Pcrturbativc amplitudes on the node fl of LAGEOS due to ocean tides. P indicates the periods in 
days, A the ampUtudes in mas and E the percent error in the C;^j- The values employed for them and the 
related errors are those quoted in EGM96 model. 



Oeean tidal perturbations on the node f2 of LAGEOS II 
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Table 7: Perturbative amplitudes on the node of LAGEOS II due to ocean tides. P indicates the periods 
in days, A the amplitudes in mas and E the percent error in the Cj^j- The values employed for them and the 
related errors are those quoted in EGM96 model. 



Oeean tidal perturbations on the perigee uj of LAGEOS II 
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Table 8: Perturbative amplitudes on the perigee ui of LAGEOS II due to ocean tides. P indicates the periods 
in days, A the amplitudes in mas and E the percent error in the Cy^^. The values employed for them and the 
related errors are those quoted in EGM96 model. 
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Mismodclcd ocean tidal perturbations on nodes of LAGEOS and LAGEOS II and perigee of LAGEOS II 
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5S1^ -* /Ari"* "* (%) 


Suj^^ (mas) 




Sa i=2 p=l 9=0 


6.7 


1.37 


1.1 


2.5 


1.9 






Sa 1=3 p=l 9=-l 


10 










11.4 


5 


Sa 1=3 p=2 9=1 


10 










29.7 


13 


Ssa 1=3 p = l 9=-l 


27.2 










6.2 


2.7 


Ssa 1=3 p=2 9=1 


27.2 










9.8 


4.3 


Ki 1=2 p=l 9=0 


3.8 


5.9 


4.7 


1.3 


1 


6.76 


2.9 


Ki 1=3 p=l 9=-l 


5.2 










64.5 


28.3 


Ki 1=3 p=2 9=1 


5.2 










18 


7.9 


Ki 1=4 p=2 9=0 


3.9 






1.6 


1.2 






Pi 1=3 p=l 9=-l 


18.5 










5.3 


2.3 


Pi 1=3 p=2 9=1 


18.5 










6.4 


2.8 


K2 1=3 p=l 9=-l 


5.5 










11.7 


5 


K2 1=3 p=2 9=1 


5.5 










4.8 


2 


S2 1=3 p = l 9=-l 


7.1 










6.9 


3 


S2 1 = 3 p=2 9=1 


7.1 










4.4 


1.9 


T2 (=3 p=l g=-l 


50 










2.6 


1.1 



Table 9: Mismodeled ocean tidal perturbations on nodes O of LAGEOS and LAGEOS II and the perigee u) 

of LAGEOS II compared to their gravitomagnetic precessions over 4 years. The effect of the ocean loading 
has been neglected. When the 1% cutoff has not been reached a - has been inserted. The values quoted for 
Ki I = 3 p=l includes also the mismodeling in the ocean loading coefficient ^3 assumed equal to 2%. 



Solid Earth even zonal tidal contributions to S/i 


Tide 


1=2, m=0, ] 


3=1, q=0 


A{^lj) (mas) 


A{Qji) (mas) 


A{uJii) (mas) 


61J, 


055.565 


-1079.38 


1982.16 


-1375.58 


-0.219 


055.575 


5.23 


-9.61 


6.66 


1.06 • 10-3 


056.554 Sa 


9.95 


-18.28 


12.69 


2.02 • 10"3 


057.555 Ssa 


31.21 


-57.31 


39.77 


6.33 • 10-3 


065.455 


5.28 


-9.71 


6.74 


1.07-10-3 


075.555 Mf 


4.94 


-9.08 


6.3 


1 • 10-3 



Table 10: Contribution of the even zonal solid tidal constituents to Sfi by means of the formula 5fl^ + SCl^^ x 
0.295 - 6lu" X 0.35 = 60.05 x /x. 



Ocean even zonal tidal contributions to 


Tide 


1=2, m=0, p=l, q=0 


A{Qi) (mas) 


A(Qii) (mas) 


A(oi)ii) (mas) 




056.554 Sa 


-20.55 


37.71 


-26.17 


-3.68-10-3 


057.555 Ssa 


-5.98 


10.98 


-7.62 


-1.28 - 10-3 


065.455 Mm 


-0.54 


1 


-0.69 


-8.78 - 10-5 


075.555 Mf 


-0.62 


1.13 


-0.78 


-1.73-10-^ 



Tabic 11: Contribution of the even zonal ocean tidal constituents to dji by means of the formula SCt^ +5(1^^ x 
0.295 - (5w" X 0.35 = 60.05 x /x. 
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Ocean odd zonal tidal contributions to Six 


Tide 


1=3, m=0, p=l, q=-l 


Aiyti) (mas) 


AiVLii) (mas) 


A{ujij) (mas) 


5ii 


056.554 Sa 


-0.063 


0.13 


-114.35 


0.66 


057.555 Ssa 


-9 • 10-3 


0.028 


-22.95 


0.133 


065.455 Mm 


-4 • 10-4 


1 • 10-3 


-1.53 


-8.93 • 10-3 


075.555 Mf 


-1 • 10-4 


7 • 10-4 


-0.58 


3.41 • 10-3 



Table 12: Contribution of the odd zonal ocean tidal constituents to 5^ by means of the formula 50/ + 5tl^^ x 
0.295 - Sixi^^ X 0.35 = 60.05 x /x for p = 1, q= -1. 



Ocean odd zonal tidal contributions to 


Tide 


1=3, m=0, p=2, q=l 


A{Q,i) (mas) 


A{Q,n) (mas) 


A{ijJn) (mas) 




056.554 Sa 


0.047 


-0.36 


297.34 


-1.72 


057.555 Ssa 


7 • 10-3 


-0.044 


36.07 


-0.209 


065.455 Mm 


4- 10-4 


-2 ■ 10-3 


1.642 


-9.55 ■ 10-3 


075.555 Mf 


1.79 • 10-4 


-7-10-4 


-0.60 


-3.52-10-3 



Table 13: Contribution of the odd zonal ocean tidal constituents to by means of the formula 50,^ + SCl^^ x 
0.295 - (5w" X 0.35 = 60.05 x /x for p = 2, q=l. 



